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Abstract 

In light of the recent swampland conjectures, we explore quantum cosmology and eternal inflation 
beyond the slow roll regime. We consider a model of a closed universe with a scalar field ¢ in the 
framework of tunneling approach to quantum cosmology. The scalar field potential is assumed to 
have a maximum at ¢ = 0 and can be approximated in its vicinity as V(¢) + 3H? — 5m". Using 
the instanton method, we find that for m < 2H the dominant nucleation channel for the universe 
is tunneling to a homogeneous, spherical de Sitter space. For larger values of m/H, the most 
probable tunneling is to an inhomogeneous closed universe with a domain wall wrapped around its 
equator. We determine the quantum state of the field ¢ in the nucleated universe by solving the 
Wheeler-DeWitt equation with tunneling boundary conditions. Our results agree with earlier work 
which assumed a slow-roll regime m < H. We finally show that spherical universes nucleating with 
m < 2H undergo stochastic eternal inflation with inflating regions forming a fractal of dimension 
d > 2. For larger values of m the field ¢ is unstable with respect to formation of domain walls and 


cannot be described by a perturbative stochastic approach. 
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It is a widely held view that in order to understand how the universe originated we should 


treat the universe quantum mechanically and describe it by a wave function V rather than by 


a Classical spacetime. The picture that has emerged from this approach is that a small closed 


universe can spontaneously nucleate out of nothing, where ‘nothing’ refers to a state with 


no classical space and time[1—6]. The cosmological wave function V can be used to calculate 


the probability distribution for the initial states in the ensemble of nucleating universes. 


At least some of these universes are expected to go through a period of inflation, driven 


by the potential energy V(@) of the inflaton field ¢. This energy is eventually thermalized, 


inflation ends, and from then on the universe follows the standard hot cosmological scenario. 


Inflation is a necessary ingredient in this kind of scheme, since it provides the only way to 
get from a tiny nucleated universe to the large universe we live in today. 


The wave function of the universe W satisfies the Wheeler-DeWitt (WDW) equation 
1Y =, (1.1) 


where H is the Hamiltonian operator. To solve this equation, one has to choose some 
boundary conditions for UV. The most developed proposals for this choice are the Hartle- 
Hawking [2] and the tunneling [7, 8] boundary conditions. The predictions resulting from 


the two proposals have mostly been studied assuming a slowly varying potential V(d),' 
W/V<«l, |V"\/V«l (1.2) 


in reduced Planck units (87G = 1). Then the Hartle-Hawking wave function predicts that 
the universe is most likely to nucleate with the lowest (positive) value of V(#) and thus 
disfavors inflation. The tunneling wave function, on the other hand, predicts nucleation with 
@ near the maximum of V(@), which is a favorable initial condition for (hilltop) inflation. 
Here we shall adopt the tunneling proposal. 

With this choice, and assuming that the conditions (1.2) are satisfied, most of the uni- 
verses nucleate with ¢ close to the maximum of V(@) and immediately enter the regime of 
stochastic eternal inflation [9-11]. The dynamics of the field ¢ near the top of the potential 
is dominated by quantum fluctuations and is well approximated by a random walk. In any 
horizon region the field gradually drifts away from the top and eventually rolls down to the 
minimum of V(@), but in different regions this happens at different times, and at any time 
there are parts of the universe that are still inflating with @ near the top. The inflating 
regions form a self-similar fractal [12], with their total volume growing exponentially with 
time. 

This scenario, however, is now being seriously questioned. There is growing evidence that 
a wide class of quantum field theories do not admit a UV completion within the theory of 
quantum gravity, even though they look perfectly consistent otherwise. Such theories are 
said to belong to the swampland, as opposed to the landscape. A number of different criteria 


that the landscape potential V(¢) should satisfy have recently been suggested [13, 14]. One 


' We assume for simplicity that the potential depends on a single scalar field ¢. An extension to multiple 
fields is briefly discussed in Section 2. 


of them is the so-called refined swampland conjecture [15, 16], which asserts that at any 


value of ¢ with V(¢) > 0 one of the two conditions has to be satisfied: 
|V'|>cV, or V"<—-&éVv, (1.3) 


where c and € are positive constants of order one. This conjecture is clearly incompatible 
with the conditions (1.2). It is also in considerable tension with the inflationary scenario, 
which usually requires a slowly varying potential V(@). 

The swampland conjecture (1.3) as it stands may or may not be true. We take an 
agnostic attitude to this issue, but the above discussion motivates us to reconsider quantum 
cosmology and eternal inflation in the regime where the slow variation conditions (1.2) are 
not satisfied. This is our goal in the present paper. We start with quantum cosmology in 
Section 2. The tunneling wave function is expected to be peaked near the maximum of the 


potential, where V(#) can be approximated as 
1 ae 
V(¢) & Vo— 5m or. (1.4) 


We have |V"|/V(0) = m?/3H?, where H = (V)/3)'/? is the expansion rate of a de Sitter 
universe of vacuum energy density Vo. The slow variation condition (1.2) is then m < H. 
Here we are mostly interested in the regime m = H. 

We start in Section 2 by discussing nucleation of the universe using the instanton method. 
For sufficiently small values of the mass parameter, py = m/H < 2, we find that nucleation 
is described by the homogeneous de Sitter instanton having the geometry of a 4-sphere. 
The initial state is then a spherical universe with ¢ = 0, as usually assumed. On the other 
hand, for w > 2 the dominant instanton describes nucleation of a ”domain wall universe.” 
In this case the initial state still has the topology of a 3-sphere, but is inhomogeneous, 
having a domain wall centered on a 2-sphere (the ‘equator’) where ¢ = 0. We briefly discuss 
possible scenarios of subsequent evolution of a domain wall universe, referring to appropriate 
literature for more detail. 

Focusing on the spherical initial configuration, we then discuss the quantum state of 
the nucleated universe in Sec.3. We expand the scalar field @(x,t) in spherical modes on 
a 3-sphere and solve the WDW equation for the wave function V(a,@,), where a is the 
scale factor and ¢, are the mode amplitudes. Of special interest is the homogeneous mode 


n= 1. We find that the probability distribution for ¢, is fully consistent with the instanton 
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analysis. For small values of ys the wave function is peaked at @; = 0, while for pp > 2 this 
behavior is modified and the top of the potential becomes a local minimum of ||. We argue 
that this marks a transition to a domain wall universe. For inhomogeneous modes, n > 1, 
the quantum state is an analytic continuation of the standard Bunch-Davies vacuum to a 
tachyonic field (with a possible exception of the dipole mode n = 2). 

We finally discuss under what conditions nucleation of the universe is followed by stochas- 
tic inflation. One expects that there is some value [tq above which stochastic inflation 
is not possible. Several attempts have been made to determine this threshold value in the 
literature with different results [17-19], so the question is still open at this time. We address 
this issue in Sec.4, with the conclusion that fimaz = 2 — that is, stochastic inflation occurs 
for all values of 4 that allow nucleation of a spherical universe. 


Our conclusions are summarized and discussed in Sec.5. 


2. INSTANTON ANALYSIS 


We consider a model of a scalar field ¢ minimally coupled to gravity. The action of the 
model is 
s= | v=Gd'e (F ~ 5(V9) - v(0)) | (2.1) 
We assume that the potential V(¢) has a maximum at ¢ = 0, where it can be approximated 
as in Eq.(1.4): 

V(¢) =3H (1 = 51) | (2.2) 
where pp = m/H. 

In the leading order of semiclassical approximation, nucleation of the universe can be de- 
scribed using the instanton method. According to the tunneling prescription, the nucleation 
probability is then given by 

Prue «x e798, (2.3) 
where Sp is the action of the Euclidean instanton solution. The dominant instanton is the 
one with the smallest magnitude of the action |Sz|. It is usually assumed [1-5] that this is 
the homogeneous de Sitter instanton with ¢@ = 0 and the geometry of a 4-sphere of radius 


H~'. The corresponding Euclidean action is 


8 2 
CS eee (2.4) 


We will now argue that this assumption is valid only for sufficiently small values of the 
mass parameter: 44 < 2. For larger values of 4 the dominant instanton is an inhomogeneous 
domain wall instanton. 

Domain wall instantons were first discussed in Refs.[20, 21] in the context of domain 
wall nucleation in the inflationary universe. They generally have O(4) x Z. symmetry and 
the topology of a 4-sphere. Here we will mostly be interested in the regime of pp *& 2, 
corresponding to the transition between de Sitter and domain wall instantons. In this 
case the instanton solution covers only a small range of ¢ near ¢ = 0, so we can use the 
approximate form of the potential (2.2) and the geometry is accurately described by the 


Euclideanized de Sitter metric 
ds? = H~*{dC? + cos? C[dx? + sin? x(d6? + sin? 6d¢”)]}. (2.5) 


The domain wall instanton can be found with the ansatz ¢@ = ¢(¢). The field equation for 
@ is then 

o” —3tan¢¢' + wo =0 (2.6) 
with boundary conditions ¢/(—7/2) = @'(a/2) = 0. With uw & 2, this equation has two 
solutions: ¢ = 0 and 


o(¢) » AsinGg (2.7) 


with A = const. The first solution corresponds to de Sitter instanton and the second to 
domain wall instanton. 

The constant A depends on the nonlinear terms in the potential V(@) [21]. Its particular 
value will not be important for our discussion, but as an illustrative example let us include 


in the potential (2.2) an additional quartic term: 
i 
V(¢) = 3H? — see + r¢*. (2.8) 


Substituting this and the solution (2.7) into the Euclidean action 


= 1 
Se= 2H | d¢ cos? ¢ 50? + vio) ; (2.9) 
0 
we find 
Qn? 2 4 
= (3H =F = 4) A? | | 2,10 
Se = pa [BH — PW? ~4)A? + AAS (2.10) 


This action is minimized by 
7 
eee eee 2. 
ax) (2.11) 


for yp > 2 and by A = 0 for yw < 2. Note that the domain wall instanton solution (2.7), 
(2.11) does not exist for js < 2. Note also that this threshold value depends only on jz and 
is not sensitive to the nonlinear part of the potential. A numerical calculation in Ref.[21] 
confirms that seen) < se) for all 2 > 27 

Negative modes of the de Sitter instanton have been discussed in the Appendix of Ref.[21]. 
It shows that a homogeneous negative mode ¢ = const is always present for pz? > 0, as it 
should be in a tunneling instanton. There are no other negative modes for  < 2, but for 
jt > 2 additional, dipole negative modes of the form (2.7) appear. They signal that the de 
Sitter instanton is no longer dominant and that the leading instanton has a domain wall 
structure. 

We conclude that the dominant nucleation channel in models with yp > 2 is to a domain 
wall universe. Possible scenarios for subsequent evolution of such a universe depend on the 
shape of the potential V(¢) away from ¢ = 0; they have been discussed in Ref.[25]. A few 
characteristic examples are illustrated in Fig.1. 

In Fig.la the potential is a double well with a maximum at ¢ = 0 sandwiched between 
two symmetric de Sitter minima. Then at late times we have an approximately spherical 
expanding de Sitter universe with a domain wall on its equator. An important point is that 
domain walls of a fixed physical thickness exist in de Sitter space only for pz > 2 [26, 27]. 
In the opposite case, wall-like configurations like Eq.(2.7) are smeared by the expansion of 
the universe, evolving towards a homogeneous ¢@ = 0 solution. In Fig.1b the minima of the 
potential are at V = 0. The universe then evolves to an inflating domain wall solution with 
two asymptotically Minkowski regions bounded by an expanding domain wall [28]. More 
precisely, the regions on the two sides of the wall have the geometry of open (k = —1) 
FRW universes. The minima of V(¢) do not have to be symmetric; an example is shown in 
Fig.1c. Once again, the wall separates two regions containing open FRW universes evolving 


with @ > 0 and @ < 0 branches of the potential. With a suitable choice of parameters, the 


? This is similar to Coleman-de Luccia instanton [22] describing bubble nucleation in false vacuum. This 
instanton solution is dominant when it exists, but it disappears and is replaced by the Hawking-Moss 
instanton [23] when the tunneling barrier becomes sufficiently flat near the top [24]. 


co) 


= 4 0 1 7 ¢ al 
(a) Two de-Sitter vacua symmetric about (b) Two Minkowski vacua symmetric about 


p= 0. @= 0. 


(c) An anti de-Sitter vacuum at ¢ < 0 anda 


Minkowski vacuum at ¢ > 0. 


FIG. 1: A plot of V(@) with w > 2 and a local maximum at ¢ = 0 in between two local 


minima. 


potential can have a flat portion which can lead to a long inflationary period; an example 
is discussed in [25]. The potential minimum at ¢ < 0 in Fig.1c is at V < 0. As a result the 
FRW universe on that side of the wall develops a big crunch singularity [29]. 

We finally briefly discuss multi-field models where the potential can be approximated 


near its maxima as a 
1 
_ ss 242 
Vio) =3f . ze m2 ;. (2.12) 


Instanton solutions for such models in the de Sitter background have been discussed in 
[21] assuming that the potential has an O(V) symmetry, so that all mass parameters m? 


are equal to one another. For N = 2 the de Sitter instanton dominates if m; < 2H. 


Otherwise, the dominant nucleation channel is to a ‘global string universe’, having a global 
string wrapped around its equator. Similarly, for N = 3 with m; > 2H, the dominant 
tunneling is to a ‘global monopole universe’, having a global monopole and anti-monopole 
at its opposite poles. For models with lower symmetry the nucleated universe would have a 
more complicated defect structure. For example, in an N = 2 model with Z2 x Z. symmetry 


we would have m; # mz and the global string would be attached to four walls. 


3. QUANTUM STATE OF THE UNIVERSE 
3.1. General formalism 


We shall first review the standard minisuperspace analysis of the tunneling wave function. 
It was originally performed in the slow-roll regime pw < 1 (e.g., [30, 31]). Here we will extend 
this analysis, but we will keep ?¢? < 1, so that the expansion about the maximum of the 
potential is justified. 


The metric is assumed to be homogeneous, isotropic and compact: 
ds” = a(n)? (—dn? + dQ3) (3.1) 


where a(7) is the scale factor, 7 is the conformal time and dQ3 is the metric on a unit 
3-sphere. 
We expand our field @ in spherical harmonics Q: 


b= (20)? $7 brim (t) Qin (X)s (3.2) 


n,l,m 


where (27?)!/? is a normalization factor and n, 1, m are harmonic numbers on a 3-sphere. 
For further convenience we will denote all harmonic numbers as n. 
The Wheeler-DeWitt (WDW) equation for the model is obtained through the quantiza- 


tion of the Hamiltonian constraint. In this case it recasts to 


1 @& 2 
542 da 6n°V (a) — dH U(a, on) = 9, (3.3) 


where we have neglected the factor ordering ambiguity and 


V(a) =a’ (1— H’a’), (3.4) 


te ME ws 
n= 2a2 a¢? a Fawn Pn: (3.5) 


w? = (n? — 1)a* — m7a°, (3.6) 


where n € {1, 2, 3...}. 
The next step is to solve the WDW equation. We will treat the scalar field modes @, as 


small perturbations. Then a solution to (3.3) can be obtained with the ansantz 


n 


WV x exp 120259) — 5Ra(a)o? ; (3.7) 


Substituting the above to (3.3) and neglecting subdominant terms in the WKB expansion 


along with terms of order O (4), we have 


(2) ~V(a) =0 (3.8) 


dS, dR. 
2 0 n\_ p2 2 ye 
a (=) ( P ) Ri t+w, = 0. (3.9) 


The solution to (3.8) is straightforward and given by 


and 


dase is 1 — Hq?)°? 
ae = +/V(a) — S¥(a) = +! TE ) (3.10) 


where 55 (a) correspond respectively to the decreasing and growing branches in the under- 
barrier region Ha < 1. In the classically allowed region Ha > 1 they are the contracting 
and expanding branches. 


We now turn our attention to (3.9). We make a change of variables 


Via <0) 
“ Je) (3.11) 
v —,/V(a) T>0 
where 7 is the Euclidean conformal time related to 7 through 7 = 77 and 
1 
= ———__.. 3.12 
“A cosh(T) ea?) 
With this change (3.9) becomes 
dR= 12 
A RE + we = 0. l 
ie Rr +w7 =0 (3.13) 
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This is a Riccati equation and can be linearized via the transformation 


a? dy 
Ri(r) = -——. 3.14 
n (T) ce (3.14) 
The resulting ODE recasts to 
iy + 2-0 — (n? —1— ma?) m = 0, (3.15) 
a 


where we substituted the explicit form of wy. 


The solutions to the above can be expressed in terms of the associated Legendre functions: 
V»(T) = A, coshT P',(tanhr) + B, cosht Q*,,(tanhr), (3.16) 


where >» oe 2 

Y= € + mm) 5 (3.17) 
and A,, B,, are constant coefficients. The functions v,, correspond to the mode functions of a 
tachyonic field db of mass im in the de-Sitter background geometry. Selecting the integration 
constants in (3.16) is equivalent to selecting the quantum state of the scalar field d. 

The tunneling boundary conditions yield a total of three branches: growing and decaying 
under the barrier and expanding in the allowed region. The under the barrier components V+ 
and W_ are to be evaluated for 7 < 0 and 7 > 0 respectively. The expanding branch in the 
Lorentzian region can be found by analytically continuing V_ to imaginary conformal time 
n = it. The coefficients multiplying the exponentials in Eq.(3.7) can then be determined by 
matching at the turning point T = 7 = 0 [30]. 


3.2. Regularity conditions 


In order to specify the tunneling wave function we will need to determine the arbitrary 
coefficients A,, B, in Eq.(3.16). This can be done by requiring that the wave functions 
for the mode amplitudes @, in (3.7) for all three branches are Gaussian rather than inverse 


Gaussian: 


R*(a)>0, Re{R,(a)} > 0. (3.18) 


nr 


Here, R* and R> are for the decaying and growing branches respectively and R,, is for 


the expanding branch. (Note that R= are real.) Following earlier literature, we will call 
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Eq.(3.18) the regularity condition. It ensures that the wave function is peaked at ¢, = 0, 
so that field fluctuations are suppressed. 
We begin by examining the behavior of the growing branch well inside the Euclidean 


region, Ha < 1.3 In the limit tT — +00 the mode functions v,, behave as* 


Un(T) Ay eO-™) + By et UT (3.19) 


where we have absorbed numerical factors in A, and B,,. For B, # 0, the dominant term is 


the second one, so that v, « e+. This would give 
dln vy, 
R-(a) = —a? — ~~ —(n+1)a? <0. (3.20) 
T 


The regularity condition can be satisfied if we set B,, = 0; then 
(a) = (n—1)a? > 0. (3.21) 


This suggests that we set 
V»(T) = cosht PY, (tanh). (3.22) 


We next check regularity at the turning point (a = H~!, 7 = 0). Continuity of the wave 

function at this point requires that 

Jim R,(7) = lim Ry(r) = Jim, Ra(n) (3.23) 
The Legendre functions and their derivatives are continuous at the origin. This implies that 
the constants A, and B,, ought to be identical and the values of R,,(7 = 0) should be the 
same for all three branches. 

It is shown in Ref.[32] that with the mode functions (3.22) the regularity condition at the 
turning point is satisfied for all inhomogeneous modes (n > 1) with pw < 2. Furthermore, it 
is shown in Ref.[33] that the regularity condition is satisfied in the entire classically allowed 
range Ha > 1 if it is satisfied at the turning point. However, the homogeneous mode n = 1 
has Ri(a = H~') < 0 for all yw? > 0. At smaller values of a it also develops a singularity 
on the dominant branch: R{(a — a,) 4 —oo at some point a, in the range 0 < a, < H7!. 


This mode therefore requires a special treatment. 


3 Here we follow the analysis in Ref.[32]. 
4 This expansion along with additional properties of Legendre functions which will be used throughout the 
paper can be found in DLMF or obtained with Mathematica. 


12 


To find the mode function suitable for the homogeneous mode, we shall use the results 
of Refs.[7, 31] for slow-roll potentials as a guide. They suggest that in the regime Ha < 1 
and pd, < 1, the relevant part of the decaying branch should behave as 
nT? Ha4 


4 


WV, (Ha <1) « exp |—6r?a? — woe, (3.24) 


where we kept the leading terms in a and ¢,. Comparing this with (3.7) we have 


2 2 FT? 4 
St(Ha <1) ® S . Ri(Ha<«1)x" 5 on, (3.25) 


where the second expression is valid for u < 1. 
On the other hand, the homogeneous mode function at tT + —oo (or Ha < 1) will 


generally behave as 
v(t + —o0) « (2A; sinyr + Byrn cosym) e7?" = Cie", (3.26) 


where we defined 
C, = 2A; sinyr + Bir cos yr. (3.27) 


For C 4 0 this gives 
odin 


dt 


which disagrees with (3.25). An agreement is achieved only if we set 


Rf (Ha <« 1) =-—a we 2a”, (3.28) 


Z 
B, = —A,—tan(yz), (3.29) 
1 


so that C; = 0. 

The value of R} (0) for the above choice of the homogeneous mode is plotted as a function 
of 4 in Fig.2. We see that regularity is satisfied for  < 2 and violated in some range above 
pf = 2. This violation of regularity first occurs precisely at fp = 2, when the de Sitter 
instanton becomes subdominant. At still larger values of 41, R{ (0) becomes positive in some 
subregions. However, inhomogeneous modes with n > 1 follow regularity patterns that 
interfere with the pockets of regularity of the homogeneous mode. This becomes clear from 
Fig.3 where we plot R,(0) for a number of leading inhomogeneous modes. When the mode 
n = 1 becomes regular at > 2, higher order modes violate regularity. 

So far we have only verified the regularity condition at a + 0 and a > H~!. By sampling 


different values of 4 numerically, we convinced ourselves that regularity holds at u~ < 2 on 
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R,(r = 0) 


20+ 


[L. 


FIG. 2: A plot of Ri(7 = 0) with respect to yw. For uw < 2 regularity holds. As larger values 


of p it still holds for some subregions. 


\ 


FIG. 3: A plot of R,(7 = 0) vs. pw for the modes n = 2,3,4,5. Higher order modes violate 


R,(7 = 0) 


30, 


20+ 


-20 | 


39 L 


regularity at larger values of yz. This leads to a pattern of interference in the pockets of 


regularity after > 2. 


both growing and decaying branches in the entire classically forbidden range 0 < a < H7!. 
There are however two exceptions to this rule. 

The first exception is the homogeneous mode n = 1 on the growing branch, which has 
R; (a) < 0 in a range of a adjacent to a = 0 for all py? > 0. The other exception is the 


dipole mode n = 2, which violates regularity on the decaying branch for 1.9 S p < 2. 
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These regularity violations are discussed in Appendix A, where we argue that they are not 
necessarily dangerous. 
Overall, our analysis suggests the following choice of the mode functions: 


coshr [P1,,(tanhr) — 2 tan(yr) Q1_,(tanhr)| n=1 
nee [Prat )— ztan(y7) Qy44( )] (3,30) 


cosh Tt P}',,(tanh rT) > 1 
We note that inhomogeneous modes follow the usual Bunch-Davies prescription, the only 
difference being that we are now dealing with a tachyonic field. As a result the modes with 
n > 1 will be regular even if we analytically continue to a field with a positive mass (ju? < 0). 
The same does not hold for the homogeneous mode which would yield 


w? Hat 
uw? |. 


Ri (Ha « 1) - ; 


(3.31) 


This is a clear characteristic of the tunneling wave function, for which the probability density 
is maximized at the maxima and minimized at the minima of the scalar potential.® 

On the other hand, the appropriate choice for the Hartle-Hawking state is to select the 
Bunch-Davies modes (3.22) for all n. This gives a regular wave function for ? < 0 — that is, 
for an ordinary (non-tachyonic) massive field. Note that the difference between the tunneling 
and Hartle-Hawking approaches is limited only to the behavior of the homogeneous mode, 


and thus they predict the same quantum state for scalar field perturbations. 


4. CONDITION FOR ETERNAL INFLATION 


Suppose now that the universe nucleates with ¢ near the top of the potential, where V(¢) 
is well approximated by Eq.(2.2) with w < 2, and has initial geometry close to de Sitter. 
We want to determine the conditions under which the universe will then undergo stochastic 


eternal inflation. 


° Note that this conclusion differs from that in Ref.[34] which discussed a free massive scalar field (u? < 0) 
model using the Bunch-Davies mode function for the homogeneous mode. This mode function is regular 
for a non-tachyonic field discussed in that paper, but would exhibit pathological behavior in more general 
models that would include tachyonic fields at the maxima of the potential. 
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4.1. From Fokker-Planck equation 


We shall first review the case of a small mass, 4p < 1, which has been extensively studied 
in the literature. The scalar field dynamics in this case is accurately described by the 
Fokker-Planck equation [11] 

Op H?O’p 1 86 f OV 
dt Bde | BHOd Gs 


(4.1) 


The field ¢ here should be understood as its average over a horizon volume ~ H~®, and the 
quantity p(¢,t)d@ has the meaning of the probability for ¢ to be in the range |[¢, 6 + dq] at 
time t. The first term on the right-hand side of (4.1) accounts for quantum diffusion of the 
field and the second term for the classical drift due to the potential. Comparing the two 
terms and using V(@) from Eq.(2.2) and ‘naive’ estimates like OV/0¢ ~ V/d, we find that 


quantum fluctuations dominate at 
ld] < by = H?/m. (4.2) 


The eternally inflating region (EIR) of spacetime can be defined as the region where ¢ is in 
this range. We note that for such values of @ 
m2 ¢ 


v; <<, (4.3) 


where the last inequality assumes that the potential energy density at the top of the potential 
is small compared to Planck density. It follows from Eq.(4.3) that the geometry of EIR is 
well approximated by de Sitter space of energy density Vo. 

Outside of EIR scalar field dynamics is nearly classical, so ¢ rolls down from the top of 
the potential and the geometry significantly deviates from de Sitter. However, this exterior 
region can affect the geometry of EIR only within a small distance ~ H~! from its boundary. 


The 3-volume of EIR at a late time ¢ can be estimated as 
Kibq 
Vernlt) ~ Yoown(3Ht) |” o(6.t)as (4.4) 
as 
where & >> 1, Vo = 27?H~? is the initial volume at nucleation and exp(3Ht) is the volume 
expansion factor of inflating regions. The parameter « depends on our definition of EIR; its 
precise value is not important for our discussion here, while the condition k >> 1 ensures 


that the integration range in (4.4) covers the region where quantum fluctuations of @ are 


non-negligible. 
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With initial condition p(¢,0) = 6(¢), the solution of the Fokker-Planck equation is given 


by a Gaussian distribution 


p(¢,t) = (20) "0" (t) exp[—¢?/20°(t)]. (4.5) 
with [12] a ; 
eS ane lex» (Fut) — 1 (4.6) 


In the limit of large ¢ it approaches a ¢-independent form 


p = const - exp(—p? Ht/3) (4.7) 
and Eq.(4.4) gives 
Verr(t) « exp(dHt), (4.8) 
where 
d=3-—p?/3. (4.9) 


It was shown in Ref.[{12] that the eternally inflating region in this case is a self-similar fractal 
of fractal dimension d. 

The solution (4.5), (4.6) was used by Rudelius [17] to determine the condition for eternal 
inflation (see also references in [17] for earlier attempts in this direction). Requiring that 
the inflating volume grows with time, he found from Eq.(4.8) the condition 4 < 3. We note 
however that the Fokker-Planck equation on which the above solution is based is applicable 
only for slowly varying potentials satisfying the conditions (1.2), which in our case correspond 
to uw <1. Here we need to explore the regime of ~ ~ 1, so we cannot rely on the Fokker- 


Planck approximation. 


4.2. From quantum cosmology 


Our approach will be based on quantum cosmology. The wave function of the universe 


in the classically allowed region Ha > 1 is given by Eq.(3.7) 


1 
WU x exp | —12275p(a) — 5 Rnla)on (4.10) 
with 
a? dy 
LoS ai 
@) =-2o (4.11) 


where a(7) = —1/H7n and the mode functions v,,(7) are obtained by analytic continuation 
T — in from Eq.(3.16). In this section it will be more convenient to use the proper time 
t = —H~'ln(—H7) instead of the conformal time 7. 

The mode amplitudes ¢, behave as independent Gaussian variables; hence the probability 


distribution for ¢, can be written as 
dP({bn},t) = [] dbnpn(bn), (4.12) 


where 
Pn(nst) = (20) ¥?o;,*(t) exp[—¢7,/207,(4)]. (4.13) 
From Eq.(4.10) it is clear that the variances o,,(t) are given by 


1 


o;,(t) = 2Re(R,(t)) (4.14) 


The mode variances ¢,, are calculated in Appendix B. 
The total field ¢, comprising all the modes, depends linearly on the mode amplitudes; 


hence its probability distribution is also Gaussian: 


dP(0,t) = p(¢, t)dd = (20) "?o" (Lt) exp[—4?/207 (t)] do, (4.15) 


where 
FSO eS Sono, (4.16) 


and n? is the mode degeneracy factor. 

As it stands, the sum in Eq.(4.16) is divergent at large n. This is not surprising: quantum 
fluctuations of @ on scale | grow as | is decreased, so a calculation of the expectation value of 
¢ at a point always gives a divergent result. This divergence can be regulated by smoothing 
the field operator d(x, t) over horizon-size regions (1 ~ H~'). We can do this using a 
Gaussian window function, or more conveniently by cutting off the summation over n in 
(4.16) at some n, ~ e”®. 

The probability distribution for a smoothed field is still given by the Gaussian distribution 
(4.15). To analyze the regime of eternal inflation, we will be interested in this distribution 
in the limit of t + oo. We found in Appendix B that in this limit the time dependence of 


Oy is the same for all n: 


o° x exp(2yHt), (4.17) 
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where 


9 .\i? 3 
a4 a 41 
y (G+) a (4.18) 
as before. Hence 7 
= >» n?o2 = CH’ exp(2yHt). (4.19) 
n=1 


The dimensionless constant C’ is estimated in Appendix B: 


~ 2 (t+ 
sin(yr) D (y+ 1) T (y +3) 


(4.20) 


In principle, C should depend on the cutoff n,. However, we show in Appendix B that 
large-n terms in Eq.(4.16) make a negligible contribution to the series, making C effectively 
independent of n,. This indicates that o? is not sensitive to the size of the smoothing region. 
We note that C' is divergent at y = 1 (uw = 2) when the modes become unstable. Another 
divergence is observed in the massless limit ~ — 0 (44 > 0) which corresponds to the infrared 
singularity in de-Sitter space (see [35-37]). 

We can now estimate the volume of EIR in the same way as we did in the preceding 
subsection: 


Verr ~ Vo T p(o t)d@ = Vy 2h Oa g-1Ht (4 21) 
-Kby HVC . 


The result depends on the definition of EIR (parameter «), but the fractal dimension d of 


the inflating region depends only on the time dependence of the inflating volume: 


VerrR ie eet = ee. (4.22) 
This gives 
9 9 : 1/2 
d=3-y==-|{- : 4,23 
v=5- (G+) (4.23) 


It follows that for 4 < 2 we have stochastic inflation with fractal dimension d > 2. As pu 
approaches 2 from below, the fractal dimension approaches that of a domain wall, d = 2. 
And as we discussed in Sec.2, for 4 > 2 the dominant instanton describes nucleation of a 
domain wall universe. 

We derived the fractal dimension (4.23) assuming that the universe originated by quantum 
nucleation from nothing. This result, however, is more general; it is not sensitive to the 
assumption about the initial state. For example, Guth and Pi [38] calculated the scalar field 


variance in de Sitter space with flat spatial sections assuming a thermal initial state. They 
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found the same asymptotic time dependence of (7) as in (4.19), and this would give the 
same fractal dimension of EIR. as in (4.23).© This conclusion would formally also apply for 
p> 2, suggesting d < 2. However, in this case de Sitter space admits stable domain wall 
solutions, so fluctuations of the field ¢ are unstable with respect to formation of domain 
walls. These are nonlinear structures; they inflate along the wall surfaces and keep a fixed 
thickness in the transverse direction. Thus a description in terms of stochastic inflation is 


appropriate only for p< 2. 


5. CONCLUSIONS 


In this paper we studied quantum cosmology and eternal inflation in the regime where 
the slow-roll conditions are violated. Our approach was based on the wave function of the 
universe with tunneling boundary conditions. We considered a model of a closed universe 
with a scalar field ¢, focusing on the range of @ near the top of its potential, where the 
potential can be approximated as V(¢) ~ 3H? — (1/2)m?¢?. We were mostly interested in 
the regime m < H. Our analysis shows that the physics of both nucleation of the universe 
and eternal inflation undergo a sharp transition at the critical value of m = 2H. 

We first used the instanton method in Section 2 to study nucleation of the universe from 
nothing. For small values of the steepness parameter 44 = m/H we find that the dominant 
nucleation channel is to a spherical de Sitter universe with the scalar field at the top of the 
potential, 6 = 0. As the value of increases and reaches the threshold = 2, another 
possible nucleation channel appears — to an inhomogeneous domain wall universe having a 
domain wall wrapped around its equator. This channel, when it exists, always dominates the 
nucleation probability. These conclusions are largely based on the work [21] in the context 
of domain wall nucleation in de Sitter space (which is described by the same instanton as 
nucleation of a domain wall universe). This earlier work, however, did not make a connection 
to quantum cosmology. The subsequent evolution of a domain wall universe depends on the 
nature of the vacua separated by the wall (which are generally different). We discussed 


possible scenarios for the cases of (A)dS and Minkowski vacua, following earlier discussion 


® We note that y < 2 was also suggested as the condition for topological inflation in which inflation takes 
place inside the cores of domain walls [89, 40]. For numerical simulations of the fractal structure of such 
models see [41]. 
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in Ref.[25]. 

The quantum state of a nucleated spherical universe is discussed in Section 3, where we 
extend earlier treatments which assumed a slow-roll potential. To this end, we expand the 
field in spherical harmonics and solve the Wheeler-DeWitt equation for the wave function 
of the universe V(a,¢,) with tunneling boundary conditions. The tunneling prescription 
instructs us to impose the regularity condition, ensuring that |W| decreases as the mode 
amplitudes @,, are increased. We find a regular quantum state for w < 2, while for larger 
values of js such states do not exist. These results are consistent with the instanton approach. 

A caveat in the above analysis is the behavior of the dipole mode which violates regularity 
in the narrow range 2 > pw < 1.9. This feature is discussed in Appendix A, where we argue 
that it is not necessarily problematic. We note that the origin of this behavior remains 
enigmatic for us and it should, perhaps, be further investigated. 

Eternal inflation in the nucleated universe is discussed in Section 4. The standard treat- 
ment using the Fokker-Planck equation is not applicable beyond the slow-roll regime. Our 
approach here is based on quantum cosmology and does not suffer from this limitation. 
The wave function of the universe provides Gaussian probability distributions for individual 
modes @,. We use this to find the distribution for the total field @ (comprising all the 
modes) smoothed over horizon regions, and then use it in turn to estimate the volume of 
the eternally inflating region (EIR). We find that for 1 < 2 this volume grows exponentially 
with time, indicating that eternal inflation does occur in this range of yw. The rate of volume 
growth is related to the fractal dimension d of EIR, and we find that d > 2 and that d — 2 
from above as ps — 2 from below. For yp > 2 stable domain wall solutions exist in de Sitter 
space, and fluctuations of the field @ are unstable with respect to formation of walls. A 
description in terms of stochastic inflation is not appropriate in this regime. 

As a final remark, we comment on the status of inflation in relation to swampland con- 
jectures. Inflationary cosmology has impressive observational support and at this time there 
seem to be no viable alternatives. It seems reasonable therefore to assume that the swamp- 
land criteria, whatever their final form will be, must be consistent with slow-roll inflation. 
It is possible for example that the constants c, c’ in Eq. (1.3) have somewhat smaller values, 
e.g., ~ 0.1 — 0.01 [42]. This would be compatible with slow-roll inflation, even though the 
possible form of the inflaton potential would be strongly restricted. Depending on the values 


of c and ¢, quantum random walk of the inflaton field and the associated eternal inflation 
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may or may not be excluded [43-45]. Some authors have even suggested that regardless of 
the validity of any specific form of the swampland conjecture, eternal inflation may be for- 
bidden by some fundamental principle [17, 46, 47]. As stated in the introduction we take an 
agnostic approach towards these views. However, let’s assume for a moment that stochastic 
eternal inflation is forbidden. The most probable scenario is then a domain wall universe 
whose walls will inflate indefinitely spawning new inflationary regions [25]. It would be in- 
teresting to explore the issues related to anthropic selection and to the measure problem in 


this new kind of multiverse. 
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Appendix A ANOMALOUS BEHAVIOR OF n= 1,2 MODES 


We noted at the end of Sec. 3.2 that the homogeneous mode n = 1 violates regularity on 
the growing branch and the dipole mode n = 2 violates regularity on the decaying branch. 
The latter violation occurs in an intermediate range of a for 1.9 S p < 2 (see Fig.4). We 
will argue that these violations are not necessarily dangerous, after discussing a related issue 
below. 

Our selection of the mode functions in Eq.(3.30) was partly based on the analysis of 
the growing branch of the wave function V_(a,@,). However, this branch is exponentially 
suppressed compared to the decreasing branch (at ¢, = 0), so keeping it while we neglect 
larger corrections to the WKB formula requires some justification. This issue was addressed 
in Refs.[30, 32] with the following argument. First note that the growing and decreasing 
branches have comparable magnitudes near the turning point, Ha ~ 1, so keeping the 


growing branch is justified in this range. Furthermore, it was argued in [32] that 


P,(a) = 5 a <1 (A.1) 


in the entire classically forbidden region Ha < 1. This means that the function V,(a, dn) 


decreases with ¢,, exponentially faster than V_(a,@,), and therefore V_ dominates at suf- 
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FIG. 4: A plot of Rj (a) for several values of jz. Lower curves correspond to higher values 


of 4, with the yellow curve which nearly violates regularity corresponding to 4 = 1.9. 


ficiently large ¢,. We thus have a continuous domain in the {a,@,,} superspace, ranging 
from a = H~! to a= 0, where W_ is non-negligible compared to V,. The inclusion of the 
growing branch can therefore be extended all the way to a = 0. The condition (A.1) was 
verified near the turning point in Ref.[32]. Here we used numerical methods to extend this 
analysis to the entire range Ha < 1. 

To illustrate the characteristic behavior of the ratio P,(a), we plotted it in Fig.5a in the 
Euclidean region Ha < 1 for the leading inhomogeneous modes with pz = 3/2. It is always 
less than 1, apart from the turning point Ha = 1, where P,(H~') = 1. By sampling different 
values of 41, we found that his behavior holds in the entire range of yu € (0,2) for all modes, 
with the following two exceptions. One exception is the homogeneous mode n = 1 which 
exhibits a segment adjacent to the turning point where P;(a) > 1 for all > 0 (see Fig.5b). 
The other is the dipole mode n = 2 which exhibits similar behavior for 1 Z 1.6 (see Fig.6). 

For modes violating the condition (A.1), we cannot trust their behavior on the growing 
branch in the classically forbidden region. For the homogeneous mode this means that 
violation of regularity by R; (a) may be spurious in this case. Note also that even if R7 (a) < 
0 in some range of a < H~!, this would not necessarily be a problem. The homogeneous 
mode on the growing branch behaves as the Hartle-Hawking wave function which has a 


minimum at the maximum of V(¢). This would not be problematic, as long as the growing 
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FIG. 5: A plot of P,(a) in the entire Euclidean region for the modes n = 1,2,3,4,5 and 
pe= 3/2. 


FIG. 6: A plot of P,(a) in the entire Euclidean region for the modes n = 2,3,4,5 and 
p= 1.0: 


branch remains subdominant. 

For the dipole mode, violation of (A.1) implies that we could not have used regularity 
of RZ (a) to select the mode function v2(a) for wy < 1.6. This suggests that in this range 
of yz a different dipole mode function could be selected which would satisfy the regularity 
condition. We have verified that this is indeed the case for mode functions (3.16) with a 
range of values of Bj/A 2. We note however that a specific value of Bj/Az is not selected in 
this case. The significance of the anomalous behavior of the dipole mode is not clear to us. 


This issue may require further study. 
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Appendix B FIELD VARIANCE 


In this section we will calculate the individual variances for each mode and the total 
variance of the field by determining the values of Re(R,,(a)). We will first consider the 


inhomogeneous modes n > 1. An explicit calculation yields: 


Ha? 
Re(R, = — —_____ B.1 
e(Ba(@)) = pe (B.1) 
where the functions h;, i = 1,2,3 are given respectively by: 
P(-l1-y+n) 
hy ~ (Ha)?t?7 ; B.2 
TN ea eG) Be) 
ho ~ 1, (B.3) 
P(2+y+n) 
hs ~ (Ha)~ 8+?) B.4 
where we factored out O(1) terms independent of n and a. 
For n >> 1 we can use the asymptotics of Gamma functions to find 
hy ~ | — B.5 
(Sy (B5) 
n \3+27 
hs~(z-) B.6 
s~ (4% (B.6) 


(Note that the range of interest for y is 0 < y < 1.) It follows that forn < n, ~ Ha we 
have hy >> he, h3; hence we can drop hz and hg in Eq.(B.1). Then using Eqs.(4.14),(B.5) we 
obtain 


ow Hn *l(Ha)y™ (B.7) 


Since y > 0, successive terms in the sum of Eq.(4.16) decrease faster than n~', and thus the 
sum is not sensitive to the cutoff n,. For n ~ n,, we have hy ~ hg ~ hz, but this does not 
change the estimate (B.7) by order of magnitude. In any case, the estimate of o,, at large n 
is unimportant, since the large-n terms give a negligible contribution to the sum. 
Including now O(1) factors, we find 
ie aD (y+2+n) (Ha) 
"eH —l+miPO+9r # 


(B.8) 


where the absolute value is necessary to include the homogeneous mode n = 1. To a good 


approximation we can extend the summation in Eq.(4.16) to n > oo. This gives 
oe CH? (Ha), (B.9) 
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